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Abstract 

GS ! There are well-known bounds on light neutrino masses from cosmological energy density arguments. 

Os . These arguments assume the neutrinos to be non-degenerate. We show how these bounds are affected 

if the neutrinos are degenerate. In this case, we obtain correlated bounds between neutrino mass and 
degeneracy. 



OS 



- 1—1 

X 



In the standard big bang cosmology, neutrinos play a very important role in the evolution of the universe 
jD. As a result, neutrino properties can be significantly constrained from various cosmological parameters. 
Perhaps the most famous of such constraints is the one on the mass of light neutrinos p|. Using modern 
values of the cosmological parameters some of which we will discuss later, this bound comes out to be 



^ rrn < 46 eV , (1) 
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qq \ The derivation of this bounds assumes that the cosmological neutrinos are non-degnerate, i.e., the number 

' of neutrinos and antineutrinos are equal. If this condition is not satisfied, the bound is modified. For example, 
if the neutrinos are assumed to be massless and at zero temperature, the degenerate Fermi gas of neutrinos 
can oversaturate the energy density of the universe unless the chemical potentials \ii associated with the 
^ ■ neutrinos Vj, satisfy the constraint pj 



where mi is the mass of the light neutrino 



E 



1/4 



<7.4xlO~ 3 eV. (2) 



The two cases discussed above are extreme cases — one a bound on mass for vanishing chemical potential, 
and the other a bound on the chemical potential for vanishing mass. More generally, both mass and chemical 
potential may be non-zero. One then should obtain correlated bounds on mass and degeneracy of the 
neutrinos. In view of the importance of the neutrinos in the physics of the early universe, such bounds are 
potentially important. The purpose of this paper is to derive such bounds. For the sake of simplicity, we 
assume that only one light neutrino species dominates the energy density of the universe. 
To discuss these bounds, it is useful to introduce the parameter F u , defined by 

F„ (3) 
Po 

where the numerator is the sum of the present energy densities of the neutrinos and antineutrinos, and po is 
the total energy density of the present universe. Obviously, F v < 1. Modern theories of structure formation 
in the universe prefer Q] a value around 0.25 for this parameter. We will present our results for various 
values of this parameter. 
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The total energy density po is not very well-known. It is usually parametrized in the form 

po = 10 A h 2 n Q eV/cm 3 , (4) 

where fi is the density measured in units of the critical density, the latter being dependent on the Hubble 
parameter, which in turn is parametrized as lOO/i kms -1 Mpc -1 . Combining Eqs. (^) and (^), we can write 

(p„+p P ) a =10 4 G,eV/cm\ (5) 

where 

C„ = h 2 n F„ . (6) 

In the rest of the paper, we will take different values for this parameter £ v which measures the importance 
of the neutrinos in the early universe, and will find the combinations of neutrino mass and degeneracy which 
can produce the corresponding values for (p v + p p ) . The bounds quoted in Eqs. ([!]) and (^) correspond to 
C=0.5. 

At any epoch, the energy density of neutrinos plus antineutrinos can be written as 

Pv + pp = j Wf V^+^ [M + Hp)} 



^ J dp p 2 VP 2 + m 2 [/(p) + I{p)} , 



(7) 



where f(p) and f(p) denote the distribution functions of the neutrinos and the antineutrinos at that epoch, 
p being the magnitude of the momentum 3-vector. 

In the very early universe, the neutrinos were in equilibrium with other particles, so f(p) and f(p) 
were the Fermi-Dirac distribution functions, with appropriate temperature and chemical potential. Due to 
the expansion of the universe, there reached an epoch when the neutrino interactions were not enough to 
keep them in thermal equilibrium. This is called the epoch of neutrino decoupling. At this time, let the 
temperature of the neutrinos be To and the chemical potential po- Then, at this epoch, the energy density 
of neutrinos and antineutrinos was given by Eq. (Q), with 

exp^ V^-^ j+l 

and that of the antineutrinos can be obtained by changing the sign of pp. After that time, the neutrino 
momenta did not change at all due to collisions. They only suffered the cosmological red-shift. If the 
cosmological scale parameter was an at the decoupling era and eto now, all the momenta have decreased by 
a factor 

— aD ^ 
r = — . (9) 

The quantity r can be interpreted as the ratio of the scale factors of the era when the neutrinos were last in 
thermal equilibrium, and the present era. The distribution function of the neutrinos in the present universe 
can be obtained from Eq. (||) through the relation 

Mp) = f D (p/r). (10) 

If m = p = 0, this is still a thermal distribution with a redefined temperature rTr>. When either the mass 
or the chemical potential is not zero, the distribution is not thermal. Making a simple change of variable, 
we can write 



+ P*) = 2^L dpp 2 ^r^+m^ [f D (p) + f D (p)} . (11) 



2 



m in eV 



Figure 1: The contours of equal energy density of neutrinos plus antineutrinos as a function of neutrino mass 
and the degeneracy parameter r\ v introduced in Eq. (13). From the outer to the inner curve, the values of C, v are 



h and A respectively. The tick marks on the axes are at 2, 5 and 8 of each decade. 



Since the distribution is not thermal as we have already said, it is meaningless to talk about the chemical 
potential of neutrinos in the present universe. Rather, we can talk about the difference of the number 
densities of neutrinos and antineutrinos. Using similar argument, it will be given by 

K - = / d PP 2 [fo(p) ~ Fd(p)] ■ (12) 



In Fig. 0, we have shown for a range of neutrino mass m, the values of the degeneracy parameter 



(13) 



which give certain preassigned values for the density paramter Q v . In Eq. (|K5|), n 7 is the number density 
of photons in the microwave background which has been introduced to obtain a dimensionless parameter 
for neutrino degeneracy. We will discuss the characteristics of these results later. Before that, we want to 
discuss different considerations that go into obtaining the results. 

From the previous discussion, it seems that if we want to find out the energy density for the neutrinos, 
we need to know m, fijj, and r. The plot of Fig. |l| is two dimensional, so we need to show that two 
of these variables are dependent. We take m and /id as independent. Given these two, we need to solve 
for the temperature at which the reaction rate of the neutrinos become equal to the expansion rate of the 
universe. For this part, we can use r = 1 because, until decoupling took place, the neutrinos were in thermal 
equilibrium. Once Tjj is obtained this way, we can put 



Vrb 



(14) 



for the calculation of energy and number densities in the present universe, where To is the temperature 
associated with the microwave background. For y x ^ = 1, this follows from Eq. ([!(]) and the comment 
following it, which is valid for a Bose-Einstcin distribution as well. The issue of the departure of y r h from 
the value of unity will be discussed next. 

The value of y r h should be equal to unity provided the number of photons have not increased due to 
the annihilation of other species of particles after the neutrinos decoupled. Thus, if the neutrino decoupling 
temperature is smaller than the electron mass, we should put y r h — 1. If, on the other hand, Td > m e , we 
should take into account the reheating of photons from e~-e + annihilations. This gives 0] y r h — (4/11) 1 / 3 . 
These considerations have been taken into account in the plot of Fig. [l], as we discuss below. 



The decoupling temperature is determined by using the criterion T — where T denotes the reaction 
rate for neutrinos and H is the Hubble parameter. In the early universe, H is given by T 2 /Mp, apart from 
numerical factors of order unity. On the other hand, neutrino reaction rates are given by T = na, where 
n is the number density of the particles that neutrinos interact with, and a is the scattering cross section. 
When neutrinos are highly degenerate, they are more abundant than electrons and positrons. Therefore, 
elastic scattering as well as v-v annihilation are the most efficient mechanisms for redistribution of energy 
and momentum of neutrinos, which keep them in thermal equilibrium. Hence, in the expression for the 
reaction rate, one should use n — n v + n p . The cross section is given by G 2 F E 2 where E is the typical energy 
of the neutrinos. Since the masses are small in the range under consideration, we can approximate E by the 
temperature T . Therefore, the decoupling temperature is determined from the equation 

' (15) 



G F M P 



where n v and n e should be determined using the equilibrium distribution functions. 

The plots have been made for three different values for the density of the neutrinos and antineutrinos. 
First, considerations of the age of the universe dictates || that h 2 tto cannot be larger than about 0.5. This 
implies that < 0.5. We make the plots corresponding to three values of (i,, viz., 1/2, 1/4 and 1/8. If 
h 2 ilo is indeed 1/2, these values of corresponds to neutrinos and antineutrinos contributing 100%, 50% 
and 25% respectively to the energy density of the universe. 

The plots for a given density clearly shows three regions. First, there is a vertical region for vanishingly 
small degeneracy. Here, our bound agrees exactly with Eq. (|l]), as it should. In fact, the outer vertical line 
has an intercept on the horizontal axis exactly at 46 eV. Second, for small or vanishing mass, we obtain 
a horizontal branch which corresponds to the bound of Eq. (||). Here, the intercept on the vertical axis 
is 2.2 x 10 3 , which is the number obtained from Eq. (|2|). Finally, there is an intermediate region which 
represents the new results obtained in this paper. This part follows roughly a power law behavior, given by 



24C, . (16) 



The glitch encountered in this region corresponds to the shift of the value of y r h from (4/11) 1 / 3 to unity. 
Earlier, we explained how we have calculated the decoupling temperature Tjy. If To is less than the electron 
mass, we have taken to be unity. Otherwise, it is taken to be (4/11) 1 / 3 . Because of these sudden 
approximations, the shift appears as a glitch. In a more detailed calculation of the densities of different 
particles, involving the Boltzmann equation, there should be a more gradual shift from one value of y r h to 
another. 

The usefulness of these plots is the following. There are quite a few recent indications that neutrinos are 
massive. These include the solar neutrino data the atmospheric neutrino data Q and direct neutrino 
oscillation expermients |^|. Each kind of data points to a different range of neutrino mass, but none of them 
very close to the value obtained from cosmological dark matter considerations assuming zero degeneracy 
for neutrinos and taking values of h 2 flo close to 1/2. The point that we make here is that, if we introduce 
neutrino degeneracy as a free parameter, we can obtain solutions for any energy density given the value of 
neutrino mass. Such a solution would of course be acceptable if it conforms with other bounds on neutrino 
degneracy, e.g., those coming from primordial nucleosynthesis Q. One can also ask whether, given a certain 
value of mass, the required value of neutrino degeneracy is a plausible one. This can only be addressed 
within a certain scenario of generation of lepton asymmetry in the universe. This is a separate issue and 
should be taken up separately. 



Note added : After this paper was accepted for publication, we came to know about a paper by Khare and 
Deo flo|| which also found the correlated bounds on neutrino mass and degeneracy from the upper limit on 
cosmological energy density. However, these authors assumed the neutrino distribution to be thermal even 
in the present universe. In our analysis, we have taken into account neutrino decoupling and the subsequent 
deviation of the neutrino distribution from a thermal one. 
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